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Recent lattice results for the energy of gluonic excitations as a function of quark 



separation are shown to correspond to transverse relativistic flux tube vibration modes. 
For large quark separations all states appear to degenerate into a few categories which 
are predicted uniquely, given the ground state. 



Introduction 



Mesons in which the gluons are in an excited state have been discussed for some time. 
There are two main pictures that have evolved for treating these excited states. The first 
is the constituent gluon approach where the quarks and a gluon move in an MIT bag [M 
or a potential model |2|, ||. The second picture envisions the quarks to be connected by 
a string or flux tube |4], |^, |6j which has quantized transverse vibrations. In this case the 
flux tube can be thought of as a coherent gluonic state. In all of these models the resulting 
meson is analogous to the diatomic molecule where the gluonic degrees of freedom are the 
"electronic state" that can assume many levels of excitation. Each excited state yields an 
interaction energy that acts as an adiabatic potential in which the quarks or "ions" move. 
The ground state of the glue corresponds to standard meson states and the excited glue to 
"hybrid meson" states. 

Recently the excited glue states with fixed end points have been investigated in detail by 
lattice simulation [0. These calculations are done with an improved action in the quenched 
approximation for a variety of gluonic operators, and on several anisotropic lattices. It is our 
purpose here to point out that the systematics of the gluon states are extremely simple from 
the vibrating relativistic flux tube point of view. To a remarkable extent the gluon states 
group themselves into a few highly degenerate states at large quark separations, reflecting 
the well known degeneracy of the quantized two-dimensional harmonic oscillator. 

We further show that, given the ground state potential, the hybrid adiabatic potentials are 
uniquely predicted and agree well with the lattice results. Our calculation is fully relativistic 
and does not introduce arbitrary procedures as required by previous work IJ. 

Relativistic strings 

It is a common misconception that a free relativistic string must be formulated in twenty 
six dimensions in order to be consistently quantized. In fact, quantized theories of a single 
non-interacting relativistic string can be defined consistently in any spacetime dimension 
smaller than twenty six using the standard string theoretic methods. Long ago, Brower and 
Goddard and Thorn || showed that free bosonic string theories in spacetime dimensions 
D < 26 are free of ghost (negative norm) states as long as the first excited state is not 
tachyonic. Subsequently Rohrlich || found an oscillator quantization of the non-interacting 
relativistic string that is manifestly free of ghosts in any dimension, while Polyakov ]10[ 
quantized the string as a sum over random surfaces that is consistent in dimensions twenty- 
six or smaller. It is only in the context of dual models and their superstring offspring that the 
theory becomes consistent in a single (critical) dimension. This is because unphysical states 
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that may be consistently eliminated from a free string cannot be consistently eliminated 
from an interacting string theory ||TT . 



The Nambu-Goto action, with fixed end boundary conditions may be quantized consis- 
tently in D = 4 using the Gupta-Bleuler method in the temporal gauge, X°(a, r) = t = —r. 
The energy of a string of tension a and distance r between the fixed ends, 



E(r) = ^(ar) 2 + 2na(N + c) , (1) 

oo 

N = ]Tm(n++n-), (2) 

m=l 

follows from the zero mode of the Virasoro constraint JI2], [13| . The index m labels the mode 
level which is occupied by n+ phonons of positive helicity and n~ phonons of negative helicity. 
The constant c is an arbitrary normal ordering constant subject only to the constraint, 
c > —1, of the no-ghost theorem. 

In the temporal gauge, Lorentz invariance does not impose any restriction on spacetime 
dimension D II as it does in light-cone gauge and only requires that c be chosen such that 



the system has a rest frame. This formally infinite constant is often calculated [04, [L5L 16 



by summing the Casimir zero-point energies using zeta function regularization, yielding the 
value c = — ^fep. The standard BRST quantization method also yields the values D = 26 
and c = — 1. The Gupta-Bleuler method we use here yields a consistent quantum theory of 
a single string for any value of c > 0. It is well known that different methods of quantizing 
theories with constraints, such as the Nambu-Goto string, need not be equivalent and may 
differ in their energy spectra as well as in their dynamical degrees of freedom. 

The excited glue states are completely specified by the separation between the ends and 
the occupation number of each of the modes, = |r, {n+,n~}). These states can be 
labeled in molecular notation [17| by three quantities. The first is their angular momentum 



along the quark axis, 

oo 

A=£«-n m ). (3) 

m=l 

States with |A| = 0, 1, 2, . . . are denoted S, II, A, . . .. The second label is the CP value of 
the glue which appears as either a subscript g or u depending on whether CP is even or odd 
respectively. The £ states (A = 0) are labeled additionally by a superscript ± denoting their 
parity under reflection through a plane in which the axis lies. The CP value of the flux tube 
is determined by |4]] 

Vcp = (-If- (4) 

In the vibrating flux tube model the S + and S~ states are always degenerate as well. 

For N = and 1 the glue states are uniquely and U u respectively. For N = 2 the 
flux tube can be excited in the S 9 , Yl g , and A g states. For an arbitrary N > 1, states with 
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A = 0, 1, . . . , N can be excited. This degeneracy is a firm prediction of the flux tube picture 
and is expected to hold for large QQ separations. 

Testing the flux tube vibration picture 

In this section we compare the predictions of Eq. (|l]) for flux tube vibrations with the results 
of a lattice simulation of quenched QCD |7]]. The lattice energies are given relative to the 
ground state energy at a quark separation of r = 2r , where r is a hadronic scale 

distance determined |18[ from the E+ data at large r. 

The ground state (E+ ) potential from Eq. (|1|) with N = 0, together with a short distance 
QQ color singlet potential — 4a s /3r and an additive constant, gives 



Vg+ (r) = + VaV 2 + 2nac + C . (5) 



The existence of a potential for all values of r requires c non-negative. If c = 0, the "Cornell" 
potential is recovered, 

V ( r ) ^ V = _1^£ + C + a r . (6) 

The best fit of Eq. (|5|) to the lattice data occurs in fact where c = 0. Assuming this, we 
obtain for the best representation of 

a s = 0.234 , 
arl = 1M > 

Cr = -2.51, (7) 
c = 0. 

Once the above values have been fixed by the ground state lattice data, all of the excited 
states are predicted uniquely, with 



V N (r) = y/(ar) 2 + 2na{N + c) + C 



The gluonic excitation energies obtained by lattice simulation are shown in Figs. 1-3 in 
molecular notation. In Fig. 1 the A = states are plotted. The lowest A = 1 (IT ff)U ) 
and A = 2 (A g>u ) states are displayed in Figs. 2 and 3 respectively. For A ^ 0, the +A and 
—A states are degenerate. 

In the figures we also show the N = 1, 2, and 3 predictions superimposed upon the E, IT, 
and A lattice data 0. We observe general agreement at large r. The II and A states follow 
the vibrating string prediction even down to small QQ separations. The E~ states do not 
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closely resemble either the iV = 2 or 3 curves. One possible explanation is that they have 
admixtures of higher N excitations. 

Because the flux tube may be viewed as a coherent gluonic state, one might expect that 
for small quark separation it should continue smoothly to the perturbative configuration of 
one gluon and a QQ color octet. We test this assumption by modifying the prediction of 
Eq. (H) by adding to the excited potentials a color octet short distance repulsion, 

cv 

WO = ' (9) 

where a s was determined (|7|) by the short distance color singlet attraction of the ground 
state. Our results (dashed curves) indicate that this short range behavior seems to improve 
most of the predictions, if only marginally. 



Conclusions 

Relativistic vibrating string models have been widely discussed in the literature. Some au- 
thors maintain that a tachyonic ground state (c < 0) is required in fewer than 26 dimensions 



15| , |16| [19], . Others maintain that there are quantization ambiguities that allow, at least 
for a non- interacting string, consistent quantization in four dimensions [[5], ^|, |T^, ^T|. We 
advocate the later interpretation and show that the vibrating (phonon) excitations closely 
correspond to a lattice simulation of QCD. 

The success of this simple flux tube picture strongly supports the string vibration modes 
as being the correct low energy degrees of freedom for the gluonic excitations in hybrid 
mesons. The original work of Isgur and Paton |3| has the same large r limit as our result 
but is intrinsically non-relativistic. Since waves propagate with the speed of light on a flux 
tube, a relativistic model is mandated. Furthermore, the Isgur-Paton model introduces an 
arbitrary radial functional dependence to avoid unphysical behavior at r = 0. In this function 
a free parameter can be tuned so that the potential has zero slope at r = for the N = 1 
state. For larger N the intercept rises as required, but now the slope becomes progressively 
more negative. This effect is demonstrated numerically in ||. Our model predicts a unique 
potential for the excited states once the parameters are determined by the ground state. 
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Figure 1: The sigma states. The bottom line is the best fit of Eq. ([5]) to the states, which 
completely determines the other curves. The other solid lines are the predictions from the 
vibrating flux tube, while the dashed lines are the predictions of the flux tube plus an octet 
repulsive potential. The S g states correspond to even N, while the S n correspond to odd N. 
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Figure 2: The pi states. The bottom line is the best fit of Eq. (J|) to the states, which 
completely determines the other curves. The other solid lines are the predictions from the 
vibrating flux tube, while the dashed lines are the predictions of the flux tube plus an octet 
repulsive potential. 
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Figure 3: The delta states. The bottom line is the best fit of Eq. @ to the states, which 
completely determines the other curves. The other solid lines are the predictions from the 
vibrating flux tube, while the dashed lines are the predictions of the flux tube plus an octet 
repulsive potential. 
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